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ABSTRACT 

We use imaging obtained with the Hubble Space Telescope Wide Field Camera 3 to search for 
Z850 dropouts at z 7 and Jno dropouts at z ~ 9 lensed by the Bullet Cluster. In total we find 
10 Z850 dropouts in our 8.27 arcmin^ field. Using magnification maps from a combined weak and 
strong lensing mass reconstruction of the Bullet Cluster and correcting for estimated completeness 
levels, we calculate the surface density and luminosity function of our zgso dropouts as a function of 
intrinsic (accounting for magnification) magnitude. We find results consistent with published blank 
field surveys, despite using much shallower data, and demonstrate the effectiveness of cluster surveys 
in the search for z ~ 7 galaxies. 

Subject headings: gravitational lensing - galaxies: high redshift 



1. INTRODUCTION 

The properties of the first galaxies in the universe are 
of particular interest given their likely role in the reion- 
ization of the intergalactic medium (IGM). This phase 
change of the universe from highly neutral to highly ion- 
ized is believed to have begun several hundred million 
years after the Big Bang and ended < 900 million years 
after the Big Bang (corresponding to z > 6) (see e.g.. 
Fan et al. 2006; Larson et al. 2011). However, recent 
observations of z > 6 star-forming galaxies are not yet 
conclusive as to whether the observed bright end of the 
population could alone be responsible for reionization, or 
whether we need to invoke faint objects in higher num- 
bers or with larger stellar masses producing more rest- 
frame UV photons (e.g., Bouwens et al. 2008; Stark et al. 
2009; Bouwens et al. 2011). To answer this question, 
more high redshift observations that probe the faint end 
of the population are essential. Observational constraints 
on the fraction of photons from these early galaxies that 
escape into the IGM, and the dumpiness of the neutral 
Hydrogen (HI) in the IGM are also needed. 

Lyman break galaxies (LBGs) provide the largest and 
best-studied sample of high redshift galaxies (see e.g., 
Vanzella et al. 2009 for a recent study and Giavalisco 
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2002 for a review). Their spectra exhibit a sharp de- 
crease in energy at wavelengths below the Lyman limit 
(the Lyman Break) due to more energetic photons being 
absorbed within the LBG itself. For LBGs at z > 6, 
further absorption below Lya (the Lya Break) occurs 
in intervening HI clouds. Their high rest-frame UV lu- 
minosity blueward of Lya and distinctive spectral break 
enable detection of these galaxies out to high redshift 
through the "dropout" technique. This technique utilizes 
the significant drop in galaxy fiux as you move blueward 
of the break. LBGs have been detected at z ~ 5 as 
V-band dropouts, z ^ 6 as i-band dropouts, z ~ 7 as 
z-band dropouts, z ~ 8 as Y-band dropouts, and even 
z ^ 9 — 10 as J-band dropouts, (e.g., Giavalisco et al. 
2004; Henry et al. 2008; Stark et al. 2009; Castellano 
et al. 2010a; Bouwens et al. 2011). Standard selection 
criteria have been developed for detecting dropouts in 
a desired redshift range and avoiding interlopers outside 
that range (e.g., Giavalisco et al. 2004; Beckwith et al. 
2006; Bouwens et al. 2007, 2008). As an example of their 
effectiveness, Vanzella et al. (2009) using followup spec- 
troscopic observations find that only 10% of their z ~ 4, 
5 and 6 dropout-selected sources from the Great Obser- 
vatories Origins Deep Survey (GOODS) are interlopers 
(83% of which are galactic stars). 

Despite the effectiveness of the dropout technique, high 
redshift galaxy detection is greatly complicated not only 
by low apparent brightness due to high luminosity dis- 
tances but also low intrinsic brightness due to low stel- 
lar masses compared to moderate redshift (z ~ 2 — 3) 
galaxies (e.g., Bouwens et al. 2007). This results in faint 
observed magnitudes for sources at z ~ 7. The ability 
to detect high redshift galaxies can be greatly enhanced 
by the use of galaxy clusters as gravitational telescopes, 
as proposed by Soucail (1990), in that these galaxies are 
magnified and thus have brighter observed magnitudes. 
This method has been implemented successfully by var- 
ious authors to study galaxies over a wide range of red- 
shifts (e.g., Ahieri et al. 1999; Blain et al. 1999; EUis 
et al. 2001; Metcalfe et al. 2003; Richard et al. 2006, 
2008; Bouwens et al. 2009; Zheng et al. 2009; Bradley 
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et al. 2011) and is responsible for some of the highest 
redshift galaxies detected at the time of their discovery 
(e.g., Kneib et al. 2004; Bradley et al. 2008). 

The light from distant galaxies can be magnified by 
several orders of magnitude as a result of the large 
gravitational potential well of an intervening massive 
galaxy cluster. This magnification provides increased 
depth at the cost of decreased observed solid angle. By 
way of number counts, the effective slope of the Lumi- 
nosity Function (LF) being sufficiently steep {Peff = 
— (i(log$)/dlogL > 2) over the range of magnitudes 
probed results in a positive magnification bias in that 
reaching fainter magnitudes by means of lensing more 
than compensates for the fewer sources resulting from de- 
creased observational solid angle (e.g., Broadhurst et al. 
1995). Such steep effective slopes have been found at 
z - 7 all the way down to L - O.IL* (or M » ~ -17.8, 

Bouwens et al. 2011)^^. Going to even fainter magnitudes 
can still result in increased cumulative number counts 
as the loss at low luminosity is overcompensated by the 
steepness of the LF at the bright end. Thus the surface 
density of LBGs can increase significantly result 
of the magnification provided by cluster surveys yield- 
ing a much improved statistical analysis of the numerous 
fainter sources. 

In addition to the number of galaxies detected, magni- 
fication has the positive consequence of increased spatial 
resolution, thus allowing better examination of high red- 
shift galaxies on smaller physical scales. In blank field 
surveys of color-selected galaxies it is often impossible 
to confirm the nature of candidate high redshift galax- 
ies given the similarity of their colors with cold Milky 
Way stars. Also, due to photometric scatter lower red- 
shift (z 1 — 2) elliptical galaxies and moderate red- 
shift (z 2 — 3) dust-reddened star-forming galaxies 
can enter the high redshift selection window. In clus- 
ter surveys however lensed galaxies can be elongated, if 
they arc resolved, thus removing stars as contaminants. 
Clusters can also multiply-image background galaxies in 
a redshift-dcpendent way, thus removing lower redshift 
galaxy contaminants. For multiply-imaged systems, the 
positions of the images are redshift dependent. A z > 5 
multiply-imaged source can therefore be readily distin- 
guished from a multiply/singly imaged z 2 source and 
foreground (single) stars, provided accurate lensing maps 
are available. Cluster surveys are therefore better capa- 
ble than blank field surveys of differentiating high red- 
shift galaxies from interlopers. In addition, the magni- 
fication provided by clusters can better facilitate spec- 
troscopy of z > 7 galaxies, like that of Vanzella et al. 
2011 in which the first robust LBGs were confirmed at 
z > 7 (see also, Santos et al. 2004; Stark et al. 2007). 

There are various ingredients that make a specific clus- 
ter a better suited cosmic telescope for the study of such 
high-z galaxies. A more massive cluster produces larger 
areas of intermediate to high (/i > 2) magnification. A 
low redshift cosmic telescope would be desirable (as the 
angular area of high magnification is large) , however this 
results in undesirable amounts of obscuration of back- 
ground galaxies by cluster members. Hence clusters of 
intermediate (^ 0.3 — 0.5) redshift are better suited for 

For the LF parameterized following Schechter 1976, where 
13^ ff = -a + L/L*) 



these studies. Detailed cluster magnification maps are 
essential for the lensed galaxy analysis. These maps are 
best constrained in clusters with many strongly lensed 
images. In addition, lenses that are highly elliptical in 
the plane of the sky arc desirable because they have 
larger areas of high magnification. As a result they also 
tend to have more arcs and are thus easier to model. 
Given these considerations, cluster 1E0657-56 (Bullet 
Cluster) is in many ways ideal. Discovered by Tucker 
et al. (1995), it is very hot and X-ray luminous, highly 
elongated, and at an optimal redshift for these studies 
(z - 0.296). 

In this work we search for z > 7 galaxies in new deep 
infrared (IR) data of the Bullet Cluster obtained with 
the Wide Field Camera 3 (WFC3, Kimble et al. 2008) 
aboard the Hubble Space Telescope (HST). Following 
Bradac et al. (2009), we optimally combine strong lens- 
ing data (including the positions of the multiply-imaged 
source discussed below) with ground and space-based 
weak-lensing data using adaptive grid to reconstruct the 
magnification map of the Bullet Cluster. With this map 
we estimate the intrinsic brightness of the LBGs that we 
detect as dropouts as well as the actual observed solid 
angle of our observations. 

This paper is organized as follows. We describe the ob- 
servations and data products in Section 2. In Section 3 
we discuss our magnification maps of the Bullet Cluster. 
In Section 4 wc discuss our photometry. In Section 5 
wc outline the selection criteria employed to search for 
z > 7 galaxies, consider various possible sources of con- 
tamination, and present our final dropout sample. We 
estimate the completeness of our search in Section 6. In 
Section 7 we consider the surface density and luminos- 
ity function of our dropout sample. We summarize main 
results in Section 8. Where necessary, we assume cosmo- 
logical parameters consistent with 7-year Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) data: flm = 0.26, 
r^A = 0.74, h = 0.71, and = 0.8. AU the coordinates 
in this paper are given for the epoch J2000.0. All mag- 
nitudes are given in the AB system (Oke & Gunn 1983). 

2. OBSERVATIONS AND DATA REDUCTION 

We obtained HST/WFC3 imaging on 19-20 Novem- 
ber 2009 (Cycle 16, proposal 11099, PI Bradac), and 
with the Advanced Camera for Surveys (ACS, Ford et al. 
2003) on 12-13 October 2006 (Cycle 15, proposal 10863, 
PI Gonzalez) and 12 October 2004 (Cycle 13, proposal 
10200, PI Jones). The new WFC3 data consist of two 
WFC3/IR pointings in FllOW (hereafter Jno, 6530s 
per pointing) and F160W (hereafter Hieo, 7030s). The 
two pointings have an overlapping region centered on 
the cluster in which the above quoted exposure times 
are doubled, which are the numbers quoted in Gonzalez 
et al. (2010). The ACS/WFC imaging used here con- 
sists of two pointings centered on the main cluster. The 
main cluster has been imaged in F606W (hereafter Vgoe , 
2340s), F775W (hereafter ijrs, 10150s), and F850LP 
(hereafter zgso, 12700s), while the subcluster has been 
imaged in F435W(2420s), F606W (2340s), and F814W 
(7280s). For the high- redshift galaxy survey described 
below, we primarily use the main cluster's Vgoei 1775, 
Z850j J 110- and H160 data as described below. 

We use the Multidrizzle (Koekemoer et al. 2002) rou- 
tine to align and combine the images. To register the 
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images we determine the offsets among the images by 
extracting high S/N objects in the individual, distor- 
tion corrected exposures. We use SExtractor (Bertin & 
Arnouts 1996) and the IRAF routine geomap to identify 
the objects and calculate the residual shifts and rotations 
of individual exposures, which were then fed back into 
Multidrizzle. We use square as the final drizzling ker- 
nel and an output pixel scale of 0.1", smaller than the 
original pixel scale of the WFC3/IR CCD allowing us 
to exploit the dithering of the observations and improve 
the sampling of the point-spread function (PSF). Limit- 
ing magnitudes were estimated using 0.63" x 0.63" square 
apertures, which is approximately the same area on aver- 
age that SExtractor attributed to each object (obtained 
from the segmentation file) in our final zgso dropout cata- 
log. The 5cr hmiting magnitudes for Vgoe, 1775, 2350, Jiio, 
and H160 are 27.0, 26.4, 26.0, 26.4, and 26.1, respectively. 

We also utilize data from the Infrared Array Camera 
(IRAC, Fazio et al. 2004) onboard NASAs Spitzer Space 
Telescope (Spitzer) to assist in the discrimination of con- 
taminants in our high redshift LEG search. These data 
were obtained on 2004 December 17-18 and include imag- 
ing in all four IRAC bands ([3.6 /im], [4.5 /im], [5.8 /im], 
and [8 fim] ) with effective exposure times of 4 ks in each 
filter (sec Gonzalez et al. 2009 for details). 

3. MAGNIFICATION MAPS OF THE BULLET CLUSTER 

An accurate mass distribution of the Bullet Cluster is 
needed to obtain a magnification map from which the in- 
trinsic (true, accounting for magnification) magnitudes 
of our high redshift dropouts can be determined given 
their lensed (observed) magnitudes. This mass distribu- 
tion is reconstructed from the optimal use of multiply- 
imaged systems (strong lensing) and distortions of singly 
imaged background sources (weak lensing). The recon- 
struction is described in detail in Bradac et al. (2009) (for 
some other examples of reconstruction techniques using 
strong and weak lensing see Natarajan & Kneib 1996; 
Kneib et al. 2003; Smith et al. 2005; Diego et al. 2007; 
Jee et al. 2007; Limousin et al. 2007; Merten et al. 2009). 
The reconstruction does not assume a specific form of 
the underlying gravitational potential. It is performed 
on a pixelized adaptive grid, with more pixels nearer to 
the cluster and multiple images where the S/N of lensing 
is higher. Specifically, the course grid (in the outskirts 
where only weak-lensing data is available) is set to be 
18"/pixel. The grid is refined by a factor of 4 in the inner 
circular region (1.5' radius around the Brightest Clus- 
ter Galaxy), and by a factor of 16 in the regions where 
we see multiple images. This gives a final resolution of 
~ l"/pixel. For the reconstruction we use the ACS and 
WFC3 images described above as well as the data used 
in Bradac et al. (2006) and Clowe et al. (2006). We ob- 
tain a new model closely following the methodology of 
Bradac et al. (2009). In this model we use the catalog of 
strongly lensed images from Bradac et al. (2009), with 
the addition of the spectroscopically measured redshift 
of multiply-imaged system K (z=2.79; Gonzalez et al. 
2010). We also experiment with including the potential 
multiply-imaged system discussed in Section 5.4. Includ- 
ing this system results in a change in our intrinsic counts 
that is well within the error bars. The weak-lensing cata- 
logs are taken from Clowe et al. (2006) and are obtained 
from a combination of ACS and ground-based data that 



extend to a large FOV; we use the inner 9' x 9' here. 

Small-scale substructure not accounted for in the lens 
model (due to a finite number of multiply-imaged sys- 
tems and erroneous redshifts) will result in uncertainties 
in the magnification maps of our field. These magnifica- 
tion uncertainties will produce errors in our intrinsic (in 
absence of lensing) magnitudes and solid angle. Magni- 
fication error (ct^) maps from Bradac et al. (2009) con- 
servatively model the effects of this substructure as fol- 
lows. The main cluster and sub clusters were modeled as 
pseudoisothermal elliptical mass densities and the cluster 
galaxies were modeled as 30 Singular Isothermal Spheres 
(SISs). 50 mass clumps containing 8% of the total mass, 
which changes critical curve positions by 5 — 10" (more 
than the uncertainty with which multiple images are re- 
covered), were also randomly distributed throughout the 
mass map. 100 realizations of this cluster substructure 
were generated. The errors on the magnification were 
calculated as the standard deviation of the mean mag- 
nification at each pixel. While we produce a new z = 7 
magnification map, which includes the candidate multi- 
ple image of this work, we do not produce a new cr^ map. 
Instead wc use the map from Bradac et al. (2009) for a 
source at z = 6. This will be practically indistinguish- 
able from the map of a source at z ^ 7 due to the fact 
that the lensing strength (or angular diameter distance 
ratio between lens and source, and observer and source) 
at these redshifts for a lens at z 0.3 is nearly constant 
(see e.g., Bartelmann & Schneider 2001). 

Given that our map used a previous mass recon- 
struction we cannot use the ct^ values at the position of 
the dropouts. We also cannot use the average at the 
magnification of the dropout, because the actual may 
differ significantly from the average depending on posi- 
tion in the image plane. We therefore estimate at the 
position of each dropout by using the (t^ from the map 
at the pixel nearest the dropout where the mean magni- 
fication (from the 100 substructure realizations above) is 
within 10% of the magnification of the dropout. For our 
dropout sample the average distance from the dropout to 
this pixel is 10". These are the magnification errors that 
appear in Table 1 and that are also included in the error 
on the intrinsic magnitudes by adding them in quadra- 
ture with the photometric errors. 

Errors in the magnitudes (lensed and intrinsic) will 
give rise to errors in the number of dropouts in each mag- 
nitude bin of the surface density and LF of our dropout 
sample. As an example, there may be a magnitude bin 
that contains no counts but does contain a significant 
fraction of the error bars of the dropouts that fall in 
neighboring bins. We correct our surface density and 
luminosity function for this by using the sum of the frac- 
tions of the 1(7 error bars as the number of objects in 
each bin. 

4. PHOTOMETRY 

Galaxies were detected in the HST/WFC3 imaging 
data described in Section 2. We followed a procedure 
very similar to that of Stark et al. (2009) and Coe et al. 
(2006). We use SExtractor (Bertin & Arnouts 1996) 
in dual image mode with detection performed for our 
Z850 dropout search in a combined Jno + Higo image, 
thus improving the signal to noise and the measurement 
of the ccntroid, and for our Jno dropout search in an 
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H160 image. Photometry was performed in each of the 
Veoe, 1775, Z850, Jiio, Hieo filters. All magnitudes were 
first corrected for foreground galactic extinction given 
the B-V color excess in the direction of the bullet cluster 
(Schlegel et al. 1998), and the corresponding extinction 
values from Sirianni et al. (2005)^^ for Veoe, h75, and 
Z850 and from Schlegel et al. (1998) for Jno and Higo- 
This yields foreground galactic extinction values of (0.23, 
0.16, 0.12, 0.07, 0.05) for (Vgoe, i775, zsso, Jiio, Hieo), re- 
spectively. Aperture magnitudes were measured in fixed 
circular apertures of 0.8" diameter with SExtractor's 
MAG_APER and were used to measure colors. We cor- 
rect the aperture magnitudes for varying PSF size as de- 
scribed below. The final magnitudes in each filter were 
computed using a combination of aperture colors and the 
SExtractor parameter MAG_AUTO of the reddest filter 
(H160), i.e., 

m^ = MAG_AUTO (Hieo) + 

(MAG_APER (H160) - MAG.APER (x)) , (1) 

where x is Veoe, i775, 2350, and Jhq. 

There are various factors that need to be considered in 
our photometric analysis. First, intracluster light (ICL) 
occupies ^ 20 — 25% of our search area in Jno and Higo 
where it is brightest, as can be seen in the combined 
color image of Figure 1. Difficulties in determining the 
background given the presence of the ICL and in dis- 
tinguishing between ICL and galaxy flux may introduce 
additional errors and/or a bias in our measured mag- 
nitudes. We carefully examine the uncertainties intro- 
duced by ICL as described below. Second, the inability 
of MAG_AUTO to pick up the wings of the PSF can 
lead to measured magnitudes that are too faint. Fi- 
nally, the varying size of the PSF for the different wave- 
lengths probed by the filters must be accounted for in 
the aperture magnitudes. We investigate our photomet- 
ric errors and possible biases with detailed simulations 
in each band as follows. 

We first choose objects from the stellar locus and es- 
timate the full-width-half-max (FWHM) of the PSF, 
which we measure to be (0.16, 0.17, 0.17, 0.18, 0.19) 
arcsec in Veoe, i775, zsso, Jiio, and Hieo, respectively. 
This is slightly larger than the estimates of Bouwens 
et al. (2011) for example. This is expected given that we 
have fewer images and thus a larger multidrizzle drop- 
size (1.0), which is essentially the size of the top hat 
with which the PSF is convolved. We place 26th magni- 
tude^^ stars in the image convolved with a PSF generated 
with Tiny Tim (Hook et al. 2008) with the above mea- 
sured FHWM. Regarding the number of stars, we place 
a sufficient number uniformly in the image so that, given 
our completeness levels of Section 6, we detect ^ 100 
of them and measure their magnitudes. We do this 5 
times, combine these ~ 500 objects, and measure the 
mean and standard deviation of the distribution of mea- 
sured magnitude minus input magnitude. We perform 
3 sigma clipping from the median of the distribution to 
better measure the mean of the magnitude difference by 
excluding outliers, such as mock objects confused with 

Assuming an 05 V type star since these early star-forming 
galaxies should be dominated by young massive stars 

Chosen given that the average Higo magnitude of our zgso 
dropouts is 25.9 



real objects in the image. For MAG_AUTO in Higo we 
measure a mean of 0.15 and standard deviation of 0.16. 
For comparison, using a Gaussian PSF we find a mean 
of -0.02 and standard deviation of 0.18, demonstrating 
the effect of an extended PSF. We also do this for the 
aperture magnitudes used for our colors. With the Tiny 
Tim PSF we find a mean of (0.04, 0.06, 0.04, 0.19, 0.18) 
and standard deviation of (0.11, 0.11, 0.14, 0.09, 0.12) 
in Veoe, i775, zsso, Jiio, and Higo, respectively. Using a 
Gaussian PSF we find a mean of (0.00, -0.01, -0.01, -0.01, 
-0.01) and standard deviation of (0.10, 0.11, 0.13, 0.08, 
0.12) in Veoe, i775, zsso, Jiio, and Hieo, respectively. As 
expected, the bias is smaller at shorter wavelengths for 
the Tiny Tim PSF due to the smaller PSF size resulting 
in a larger fraction of the object flux within the 0.8" di- 
ameter aperture. The difference between the ACS and 
WFC3 PSFs is also manifested. 

For both MAG_AUTO in Hieo and the aperture mag- 
nitudes in all of the bands, we correct for the biases found 
with the Tiny Tim PSF, which is currently our best es- 
timate of the PSF shape. Even though the bias depends 
strongly on the PSF shape, a misestimate of 0.1-0.2 mag- 
nitudes would not significantly effect our our surface den- 
sity and luminosity function given the large bin sizes used 
(see Sec. 7). We also use the standard deviations found 
with the Tiny Tim PSF, which does not depend signifi- 
cantly on the PSF shape, as the minimum photometric 
error on the magnitudes. 

5. DROPOUT SELECTION 

In this section we discuss our selection of zgso and Jno 
dropouts through the use of color and S /N criteria (Sec- 
tion 5.1), intrinsic brightness cuts (Section 5.2), and in- 
frared {Spitzer and HAWK-I) data (Section 5.3). We 
then provide the properties of our final dropout sam- 
ple as well as potentially interesting contaminants (Sec- 
tion 5.4). Finally we address the possibility of contami- 
nation from point sources. 

5.1. Color and S/N Criteria 

We select our zgso dropouts by requiring all of the fol- 
lowing selection criteria; 

Z850 - Jiio > 1-0 + 0.4(Jiio- Hieo) (2) 

Z850 - Juo > 1-0 (3) 

^uo-i^ieo < 1.1 (4) 

S/N (i/i6o) > 3.5, S/N ( Jno) > 3.5 (5) 

5/iV(i775) <3, ^/7V(l/eo6) <3. (6) 

We will refer to Equations 2, 3, and 4 as the "color 
criteria" and Equations 5 and 6 as the "S/N criteria". 
The color criteria are adapted from those of Bouwens 
et al. (2008), which were developed for HST/NIC3 Jno 
and Hieo. To adjust them to be consistent with our filters 
we shift them by the amount by which the colors differ 
for a z = 7 template star-forming galaxy spectrum when 
measured with our filters versus those of Bouwens et al. 
(2008) (a difference of ~ 0.2 magnitudes). The template 
spectrum we use is from the isochrone synthesis code of 
Bruzual & Chariot (2003) with an instantaneous burst, 
an age of 100 Myr, metallicity oi Z = O.OO8Z0, and no 
dust (See e.g., Labbe et al. 2010; Castellano et al. 2010b). 



Figure 1. A combined color image of the Bullet Cluster (zgso = blue image, Jno = green image, and Hieo = red image) with 
the search area outlined {dashed white line, showing the tilt of the zgso pointing relative to the two WFC3 pointings) and the 
curve of maximum magnification of the lens model overlaid {cyan line). The zg5o dropout candidates are shown as green circles 
labelled with the numbers of Table 1. The zgso (Jiio) dropout contaminants are shown as red {magenta) circles labelled with 
the letters of Table 2. The Herschel source HLS13 of Rex et al. (2010) is shown as a white circle. The image is 3.3' x 3.3'. North 
is up and East is left. 



We refer to this template spectrum below as our "BC03 
starburst template" . 

The S/N cut of Equation 5 corresponds to the require- 
ment S/N(Jiio +-ffi6o) ^ 5(7. Deep blank field surveys 
typically impose S/N cuts of 2cr for the filters blue- ward 
of the dropout filter. However, given that we have less 
exposures to combine, we are likely to have more arti- 
facts due to cosmic rays that can cause artificial detec- 
tions. We thus impose 3cr S/N^"* cuts and then carefully 
inspect each object by eye in the Vgoe and 1775 images 
and throw it out if it is detected. We also remove arti- 
facts such as diffraction spikes and fake objects near the 
detector edges. 

The above zg5o dropout selection criteria select galax- 
ies in the redshift range 6 ^ z < 9 with a mean red- 
shift of 7.5 for objects with mnnjo — 26 (see Figure 5 of 
Section 6). The selection window is shown as the grey 
shaded region in the zgso-Jno versus Jno-Higo color Fig- 
ure 2. Sources that do not satisfy our dropout selection 
criteria are shown as green dots. The colors of low-mass 

Estimated as the flux divided by the flux error using the 
photometry diseussed in Seetion 4. 



M, L, and T dwarf stars are expected to lie in the region 
enclosed by orange lines. The typical colors of various 
galaxy types are shown as they evolve in redshift (for our 
BC03 starburst template, a dusty galaxy from Chary & 
Elbaz 2001, and an elliptical, spiral, and irregular galaxy 
from Coleman et al. 1980). 

We also search for Jno dropouts by requiring all of the 
following color and S/N selection criteria: 

Jno - ffieo > 1.1 (7) 
S/N (i/i6o) > 5 (8) 

5/iV(z850,i775,"^^606) <3. (9) 

Again, we also carefully inspect each object in the images 
and throw it out if it is detected in Veoe, 1775: or 2350. The 
above Jno dropout selection criteria select galaxies in the 
redshift range 8 < z < 11 with a mean redshift of 9.5 for 
objects with mHigp = 26 (see Figure 5 of Section 6). 

5.2. Intrinsic Brightness Cuts 

The 4OOOA break of a z ~ 1.4 (2.0) elliptical galaxy 
will be at the same wavelength as the Lyman Break for 
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0.0 0.5 1.0 1.5 2.0 
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Figure 2. zgso-Jiio versus Jiio-Hieo color-color diagram 
showing our zgso dropout selection window {grey shaded re- 
gion). The colors of our zgso dropouts with niHieo < 26.5, 
those used in the surface density and luminosity function, are 
shown [black squares) as well as those with mHmo > 26.5 
(black asterisks). Sources which do not satisfy our dropout 
selection criteria are shown in green. The redshift evolution of 
the colors of a Bruzual & Chariot (2003) star-forming galaxy 
with low metallicity {Z = 0.008^©), an age of 100 Myr, 
no dust, and attenuated for IGM absorption as a function 
of redshift according to Madau (1995) is shown with a blue 
line. The colors of a Chary & Elbaz (2001) dusty galaxy 
(Dusty,thin black line) and a Coleman et al. (1980) elliptical 
(E, red line), spiral (Sbc, dot-dashed purple line), and irreg- 
ular galaxy (Im, dashed cyan line) are also shown as they 
evolve in redshift. The orange Imes show region in which we 
expect T Dwarf colors to lie (Knapp et al. 2004). 

a z 7 (9) LBG. These elliptical galaxies as well as 
dust-reddened galaxies may enter the red portion of our 
selection window through photometric scatter. 

These lower redshift galaxies are on average intrinsi- 
cally brighter than the high z LBGs we are targeting. We 
thus exclude bright interlopers by rejecting those objects 
intrinsically brighter than the brightest LBG we might 
expect to see. We estimate the brightest intrinsic Higo 
apparent magnitude we might see as follows. A lOL* 
galaxy^^ at z = 6 (7.5), the lowest redshift allowed by 
our Z850 (Jiio) selection criteria, would have an apparent 
magnitude of 24.0 (24.2). We reject 2 zgso dropout can- 
didates and 1 Jiio dropout candidate that have intrinsic 
Hi6o apparent magnitudes niHieo i„t (obtained from the 
observed mHie„ and the magnification at the position of 
the object from the z = 7 magnification map) brighter 
than this cut. 

5.3. Infrared Data 

Whereas bright lower redshift elliptical and dusty 
galaxies can be rejected with S/N and intrinsic mag- 
nitude cuts, fainter objects are not as easily excluded. 
Contamination from faint elliptical and dusty galaxies in 
our Z850 dropout catalog is expected to be minimal given 

15 We use the Bouwens et al. (2011) best Rt z ^ 7 (8) LF to 
determine the L* used for zgso (Jiio) dropout intrinsic brightness 
cut. 



that our candidates do not have colors consistent with 
them (to at least one sigma, see Figure 2). We neverthe- 
less utilize our longer wavelength data to aid in distin- 
guishing these contaminants from high redshift LBGs. 
For our Z850 dropouts we consider the Higo - [3.6/xm] 
(Spitzer Chanel 1) color of these galaxies. Whereas our 
BC03 starburst template at z = 7 has Higo - [3.6/im] = 1, 
many ellipticals and dusty galaxies are significantly red- 
der than this. For example, both the elliptical and dusty 
galaxy templates used in Figure 2 begin to have Higo - 
[3.6/xm] > 1 at z = 1.3. By z = 3 the elliptical (dusty) 
template has Higo - [3.6^m] = 2.8 (2.1). The Higo - 
[3.6/xm] color for dusty galaxies can vary widely depend- 
ing on how dusty the galaxy is and its redshift. As an 
extreme example, we consider the dusty lensed galaxy 
behind the bullet cluster, which has Ay ~ 3.8 (Gonzalez 
et al. 2010). We can detect a galaxy such as this one, 
with HiQQ - [3.6/xm] = 3.75, at 3cr in [3.6/im] for sources 
as faint as mHieo = ^6. Whereas, some of these fainter 
ellipticals and dusty galaxies will be detected in [3.6/im], 
we do not expect our zgso dropouts to be, which is in 
fact the case for all of the objects in our sample. 

There are two objects that satisfy our Jno color crite- 
ria, S/N criteria, and intrinsic brightness cut, referred to 
below as objects J and L in Table 2 and Figure 4 (object 
K does not satisfy the intrinsic brightness cut). They are 
both quite bright; object J (L) has mnnjo i„t = 24.8 ± 0.2 
(24.6 ± 0.2) and is thus only 0.6 ± 0.2 (6.4 ± 0.2) mag- 
nitudes fainter than the intrinsic brightness cut. Both 
objects are clearly detected in [4.5/xm] in our Spitzer 
data. In [3.6/im] object J is clearly detected but ob- 
ject L can't unambiguously be separated from another 
nearby brighter galaxy. We do not expect our high red- 
shift LBGs to be detected in [3.6/im] (as discussed above) 
or in [4.5/im]. In addition, both of these objects are de- 
tected in the Ks-band of the High Acuity, Wide field K- 
band Imaging (HAWK-I, Kissler-Patig et al. 2008) cam- 
era on the Very Large Telescope (VLT) and have Higo - 
Ks ^ 0.5. Our BC03 starburst template would have to 
be at z ^ 11 in order to have an H-K color 0.5. It is very 
unlikely that these objects are at such high redshift given 
their intrinsic brightnesses. On the other hand, a dusty 
galaxy spectrum from Chary & Elbaz (2001) at z = 1.9 
(2.3) fits the Jiio, Higo, and Ks-band data of object J 
(L) quite well. These facts lead us to believe that these 
objects are most likely low redshift interlopers and thus 
we do not consider them high redshift candidates, and 
do not include them in the rest of our analysis. 

5.4. Final Dropout Sample 

Upon implementing the above zgso dropout selection 
criteria, wc obtain a final sample of 10 2^50 dropouts, 
with cutout images shown in Figure 3 and properties 
given in Table 1. The zgso dropouts are also shown as 
green circles in Figure 1. Only those zgso dropouts with 
niHieo < 26.5, where the completeness is greater than 
22% (See Section 6) are used in the surface density and 
luminosity function below. These 8 objects are shown as 
the black squares in Figure 2, whereas the 2 objects with 
niHieo ^ 26.0, are shown as asterisks. Upon implement- 
ing all the Jiio dropout selection criteria, we find no Jno 
dropout candidates. 

The majority of the zgso dropouts (7 out of 10) are 
found in the bottom right portion of Figure 1. This is 



7 



Table 1 

The Properties of zgso Dropouts 



Dropout 




ueci. 




( Jiio — Hie,o) 


[Z850 — ■JllO) 




mHi6a,mt ± Ai+p 


1 


104.65470 


-55.974464 


26.77±0.23 


0. 


.11±0. 


,18 


>1.3 


4.3±0.2 


in '1/1 


2. 


iU4.D00Z 


K cc OTi on 1 
-00. y 1 iaUi 


ZD. a i ±0.23 


0. 


.06±0. 


,18 


1 1 n_Ln 1 
l.lOztU.Dl 


6.0±0.0 


28.99_o.25 


3 


104.66736 


-55.968067 


24.97±0.16 


0. 


.46±0. 


,15 


1.40±0.48 


12±4 


27.67tgi? 


4 


104.66375 


-55.928802 


26.37±0.22 


0, 


.51±0. 


,19 


>1.3 


2.8±0.08 


27.461°:^^ 


5 


104.63437 


-55.978603 


25.91±0.22 


0, 


.07±0. 


,20 


1.32±0.59 


2.1±0.03 




6 


104.62446 


-55.951065 


25.85±0.19 


0, 


.02±0. 


,15 


>2.1 


10.±2 


28.37l!^:3f 


7 


104.64304 


-55.964756 


25.81±0.26 


0, 


.14±0. 


,15 


2.02±1.53 


5.2±0.5 


97 61 +0-28 


8 


104.64549 


-55.924828 


25.89±0.21 


0, 


.40±0. 


,16 


>1.9 


3.1±0.1 


97 1 9+0-22 
^' -^^-0.22 


9 


104.63254 


-55.963764 


26.00±0.16 


0, 


.46±0. 


,16 


1.94±1.28 


4.3±0.3 




10 


104.63015 


-55.970482 


26.37±0.16 


0, 


.10±0. 


,15 


1.57±0.68 


3.0±0.2 


27.57t°:r. 



Note. — The final error on niHigo i„t (^i+p) includes errors in the lensing (magnification) and photometry added in quadra- 
ture. Z85o-Jiio lower limits are the la limiting magnitudes computed in the vicinity of the dropout. 
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Figure 3. Cutouts of our zgso dropouts shown (from left 
to right) in Veoe, 1775, zsso, Jiio, Hieo, Jiio + Hieo, and a 
Z850J110H160 color image. The numeric labels correspond to 
the labels of Table 1. The cutouts are 7" x 7". To get an idea 
of physical size, for a source at z = 7 (9) and in the absence 
of lensing, 7" corresponds to 37 (32) kpc. For a uniform 
magnification of 12, the highest in our zgso dropout sample, 
7" corresponds to an intrinsic length of 11 (9.1) kpc for a 
source at z = 7 (9). 



likely due to a combination of small number statistics and 
the fact that the bottom right portion of the image has a 
slightly different magnification distribution according to 
our lens model. There also appears to be a detector arti- 
fact in the top left portion of each of the WFC3 pointings 
of ~ 250 arsec^ that is likely impeding detection in those 
regions. 

The properties of bright (mHigo int < 25) lower redshift 
objects that satisfy our 2350 (Jno) dropout color criteria 
but are rejected on account of S/N criteria, brightness 
cuts, or IR {Spitzer/Hawk-T) data are given in the top 
(bottom) of Table 2. Cutouts of these objects are also 
shown in the top (bottom) of Figure 4 and they are shown 
as red (magenta) circles in Figure 1. We also show the 
dust-reddened Herschel source HLS13 of Rex et al. (2010) 
as a white circle in Figure 1. 

Low- mass stars (e.g., L and T Dwarfs), supernovae 
and active galactic nuclei are possible contaminants that 
will appear as point sources in our images. We take ad- 
vantage of the high-resolution of WFC3/IR and the in- 
creased resolution due to lensing to determine which of 
our sources are resolved and which are not resolved and 
thus consistent with point sources. We do not attempt 
to estimate the sizes of those objects with mfjujo > 26.5 
given the uncertainty in measuring low surface brightness 
objects. Only those objects with mHieo ^ 26.5 are used 
below in the surface density and LF of Figures 7 and 9. 
All of these objects have measured FWHM greater than 
0.26" in Jiio and 0.21" in Hieo which is above the PSF 
size (0.18" in Jug and 0.19" in Higo)- Thus contamina- 
tion due to point sources is unimportant. 

Object I of Table 2 has some characteristics consistent 
with being a counter- image of object 3 of Table 1. It 
satisfies the zgso dropout color criteria. It has Jno — ^^leo 
and Z850 — Jiio colors consistent with object 3 within la 
error bars. Both have shapes extended in the direction 
expected by the lens model for a z > 6 source, as can be 
seen in Figure 1 in which the critical curve (cyan line) 
of the lens model is overlaid and the positions of the 
two objects are shown. Finally, given the position of one 
of the images, the mass model predicts a counter-image 
in the vicinity of the other (we have also verified that 
according to the lens model, none of the other dropout 
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Table 2 

A Partial Sample of zgso (top) and Jno (bottom) Dropout Contaminants 



Contaminant 


R.A. 


Decl. 


niHieo ± Ap 


(JllO — Hi6o) 


(^850 — Jlio) 








A 


104.69737 


-55.965012 


21.93±0.16 


0.72±0.15 


1.34±0.17 


>28.1 


>28.7 




B 


104.60057 


-55.957890 


21.64±0.16 


0.67±0.15 


1.51±0.17 


24.67±0.23 


26.81±0.37 




C 


104.61086 


-55.927860 


23.53±0.16 


0.62±0.15 


1.62±0.38 


25.24±0.25 


24.81±0.24 




D 


104.60263 


-55.957169 


20.30±0.16 


0.75±0.15 


1.46±0.17 


23.44±0.23 


25.19±0.23 




E 


104.64615 


-55.970993 


22.62±0.16 


0.86±0.15 


1.47±0.17 


25.37±0.25 


27.21±1.07 




F 


104.69501 


-55.975872 


23.14±0.16 


0.62±0.15 


1.25±0.25 


26.31±0.23 


25.86±0.23 




G 


104.62870 


-55.970612 


21.96±0.16 


0.62±0.15 


1.38±0.17 


24.64±0.23 


25.58±0.23 




H 


104.69444 


-55.972286 


22.49±0.16 


0.51±0.15 


1.53±0.17 


>28.1 


>28.7 




I* 


104.64444 


-55.965549 


24.18±0.16 


0.34±0.15 


1.36±0.26 


25.70±0.24 


27.23±0.73 


JllO 


J 


104.64898 


-55.971981 


23.12±0.16 


1.19±0.15 


2.51±0.48 


26.78±0.37 


>28.7 




K 


104.69546 


-55.969612 


23.10±0.16 


1.23±0.15 


1.45±0.63 


>28.1 


>28.7 




L 


104.63338 


-55.977394 


23.69±0.16 


1.62±0.15 


2.05±1.08 


28.16±1.58 


>28.7 



Note. — These objects have niHigo mt < 25 (except object I) and satisfy our dropout color criteria but are rejected from 
our sample as low redshift galaxies on account of S/N selection criteria, brightness cuts, or IR (Spiteer/Hawk-I) data. Similar 
photometric errors are due to assigning the minimum error discussed in Section 4. Objects A, B, and F lie very near a field edge 
(see Figure 1) and thus their photometry may be unreliable. Lower limits are the la limiting magnitudes quoted in Section 2. 



This object has mH,,„ . = 26.49 ±0.19. See end of Section 5.4. 

candidates would have observable counter- images) . 

Hovifcver, this object also has characteristics inconsis- 
tent with being a counter-image of object 3. Its intrinsic 
Hi6o magnitude of 26.49 ± 0.19 differs significantly from 
that of object 3. However, the error on the magnifica- 
tion of object 3 is likely underestimated given that it lies 
close (< 1") to the critical curve where the magnifica- 
tion changes rapidly with position (magnification errors 
are expected to be correct on average but may be under 
or overestimated for individual objects). Also, it is de- 
tected at 4.5(7 in 1775 and thus does not satisfy the S/N 
criterion of Equation 6. 

6. COMPLETENESS 
The completeness was estimated as function of redshift 
and observed Higo mangnitude (muieo) with simulations 
using the IRAF task artdata. For each Higo magnitude, 
niHiBo = (24, 25, 26, 27), we generated 100 mock galax- 
ies at 25 different equally spaced redshifts in the range 
5.5 < z < 10. The colors of each galaxy were determined 
from a template spectrum, where again we follow Madau 
(1995) to apply an intergalactic medium absorption cor- 
rection to the colors given the redshift. These galaxies 
were then added with a spatially uniform distribution 
to the Jiio -f-Hieo images as stars convolved with a Tiny 
Tim PSF scaled to have FWHM of (0.16, 0.17, 0.17, 0.18, 
0.19) arcsec in Veoe, 1775, zsso, Jiio, and Hieo, respec- 
tively. We then attempt to detect these galaxies follow- 
ing the same procedure as the actual dropouts, meaning 
the galaxies must both be detected as well as satisfy the 
dropout color and S/N selection criteria. For each muigo 
and z, the probability of detecting a galaxy as a dropout 
p(mHi6o:2)i completeness, is simply the ratio of the 
number of galaxies detected as dropouts to the number 
of dropout galaxies added to the images. This complete- 
ness was estimated with five SEDs from the isochronc 
synthesis code of Bruzual &: Chariot (2003) with instan- 
taneous bursts and the following age and metallicity (Z 
in units of Zq) combinations: 5 Myr and Z=0.02, 25 Myr 



and Z=0.02, 100 Myr and Z=0.02, 100 Myr and Z=0.05, 
100 Myr and Z=0.008. For our final completeness we 
use the average of the completeness estimates from each 
SED. 

The completeness function gives us an estimate of the 
redshift distribution of our dropout samples. Figure 5 
shows p (muieo = 26, z) for our zgso and Jno dropout 
samples. The completeness is also used to correct the sur- 
face density of dropouts for incompleteness. To do this 
at a specific magnitude, say mH^eo = t^o, we correct the 
surface density using p (mg), the probability with which 
we detect dropout galaxies that have the redshift distri- 
bution given by p (mg, z). We find p (mHieo ) = 58%, 54%, 
40%, and 4% for mnieo = 24, 25, 26, and 27, respectively. 
Once again this is the average of the completeness using 
each of five SEDs. As an example of how this varies with 
SED, the completeness at 26th magnitude ranged from 
37% to 45% depending on the SED. Finally, p(mHi6n,z) 
is also used to correct the LF of dropouts as described in 
Section 7. 

The ICL (occupying 20 — 25% of the search area) 
and obscuration by cluster members (obscuring ^ 25% 
of the area in the regions of the search area not near the 
ICL) are the two main factors that lower our complete- 
ness levels. Nevertheless, it is still advantageous to use 
clusters (like the Bullet Cluster) in the search for high 
redshift galaxies given that their magnifications facilitate 
reaching the more numerous faint sources overcompcn- 
sating for the loss of search area obscured by the ICL. 
The above completeness is an average estimate over the 
entire search area, including the ICL. 

We have estimated the completeness as a func- 
tion of Icnsed magnitude, p(mHi(5o), above. However, 
the completeness as a function of intrinsic magnitude, 
p(mHi60i„t)7 is required to estimate the intrinsic sur- 
face density. This is the probability that a galaxy in 
the source plane with intrinsic magnitude mfjieo i„t will 
be detected and satisfy our dropout color and S/N crite- 
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Figure 4. Cutouts of the partial sample of dropout contam- 
inants of Table 2. These objects have mnmo mt < 25 (ex- 
cept object I) and satisfy the zsso (top) and Jiio (bottom) 
dropout color cuts but do not satisfy the other criteria. They 
are shown (from left to right) in Veoe, hrs, zgso, Jiio, Hieo, 
Jiio -I- Hi6o, and a zgsoJuoHieo color image. The labels cor- 
respond to the labels of Table 2. The cutouts are 7" x 7". 

ria. It could be estimated by populating the source plane, 
properly lensing each object to the image plane, and esti- 
mating the completeness. This is difficult however given 
that the objects will be multiply-imaged. Therefore, we 
estimate p(niHi60int) the product of the probability 
p (^) that a galaxy in the source plane will be magnified 
by fi and the probability p (muieo ) that if magnified by 
/i it will be detected as a dropout. This product is then 
summed over all possible magnifications from the lens 
model. This is given by 

P (mHieoant) = '^PMPi'^liieo) : (10) 

where mniso = mHi6o,.nt - 2.51og(Ai). 
We determine p (/x) (the distribution of /x values in the 
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Redshift 

Figure 5. For an LBG with mnigo = 26, the probability of 
both detecting it and it satisfying our color and S/N criteria 
for Z850 {solid line) and Jno {dashed line) dropouts is plotted 
as a function of redshift. The probability is averaged over 
5 Bruzual & Chariot (2003) SEDs discussed in the text with 
varying ages and metallicities and then smoothed with a box- 
car average over 5 points to better represent a more diverse 
set of spectra. 

source plane) from the distribution of /i values in the im- 
age plane with each value weighted by the source plane 
area that it occupies (weighting l//i)^^. This function is 
plotted in Figure 6 as a function of magnification in mag- 
nitudes. Following Equation 10 we find p(niHi60i„t) — 
56%, 46%, 21%, and 5% for mnieoant = 26, 27, 28, and 29, 
respectively. The completeness as a function of intrinsic 
magnitude and redshift, p (mnieo int ; 2) , used to correct 
the intrinsic LF below is also estimated following Equa- 
tion 10 with p (mHiso) replaced by p (mHigo, have 
performed a check by populating the source plane and 
ray-tracing objects with a simple lens model to confirm 
that the approach of Equation 10 is accurate. 

7. SURFACE DENSITY AND LUMINOSITY FUNCTION 

In this section wc discuss our "lensed" (observed) and 
"intrinsic" (true, accounting for magnification) counts in 
terms of the surface density and LF. We also compare 
our counts with other surveys and discuss implications. 
Only those objects with lensed magnitudes ^Hieo < 26.5 
where the completeness is > 22% are considered here. 
This includes all objects in Table 1 except objects 1 and 
2. 

The completeness corrected surface density of our 
dropout sample is shown in Figure 7 for both lensed 
and intrinsic counts. We first discuss the lensed sur- 
face density. The lensed bins are 0.8 magnitudes wide 
and are centered at 25.4 and 26.2. The lensed or image 
plane area effectively observed as a result of the com- 
pleteness is given by A^ff (mjjj^gg 

60 ) ^imgi 

where Qimg = 8.27 arcmin^ is the image plane solid an- 
gle (outlined with the white dashed line in Figure 1). 
^oft.ion (111H160) applied to each bin by dividing each 

This docs not account for multiple imaging of the source plane, 
which we find through simulations to be a small cfEcct (see Sec- 
tion 7). 
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2.5 Log,o(M) 

Figure 6. The probability, according to our lens model, that 
a galaxy placed randomly in the source plane will be mag- 
nified by 2.5 log(/^) magnitudes. The function is normalized 
such that = 1- 

dropout within the bin by Acsjcn (mHieo) linearly inter- 
polated at the lensed magnitude of the dropout. 



source position. The cluster lens not only magnifies the 
brightness of the background galaxies but also magnifies 
the field of view, thus decreasing the observed solid angle. 
We estimate from the lens model a fractional decrease in 
solid angle of 0.19 ±0.02. The uncertainty was estimated 
by evaluating the solid angle change using the map of 
magnification errors described in Section 3. While using 
the magnification maps for a different LEG redshift does 
change the positions of the critical curves (by ^ 5" be- 
tween z = 5 and z ^ oo, larger than the expected accu- 
racy of the mass model), it changes the total solid angle of 
the field only negligibly (within the errors given above). 
The intrinsic or source plane area effectively observed as 
a result of the completeness (see e.g., Richard et al. 2008) 
is given by A 

cff.int (niHieo.int) — P (^Hi6o,int) ^src5 where 

Qsrc = 1-57 arcmin^ is the source plane solid angle. This 
function is shown in Figure 8. Agffjnt (mRieo int) "^^^ ap- 
plied to each bin by dividing each dropout within the 
bin by ^eff.int (niHieo int) linearly interpolated at the in- 
trinsic magnitude of the dropout. Our determination of 
the intrinsic surface density and LF does not account for 
some of the sources being multiply-imaged (i.e., a dou- 
bly imaged source is not counted twice). We have how- 
ever verified with ray tracing, that in counting only the 
brightest source in the observational setup, the expected 
number counts change by only 15%, which is well within 
the error bars as discussed next. 
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Figure 7. Lensed {triangles) and intrinsic {circles) complete- 
ness corrected surface density of zgso dropouts with bins of 
width 0.8 magnitudes. The lensed bins are centered at 25.4 
and 26.2 and the intrinsic bins are centered at 27.1, 27.9, and 
28.7. The counts are divided by the appropriate factor for 
comparison with data with 0.5 magnitude bins. The com- 
pleteness corrected z ~ 7 counts of the HUDF09 blank field 
survey of Bouwens et al. (2011) are also shown for comparison 
{X symbols). We do not however expect our intrinsic counts 
to agree with those of Bouwens et al. (2011) given that we use 
different filters and thus have different redshift distributions. 
Uncertainty in our completeness estimate (see Section 6) is 
not reflected in the error bars. 

For the intrinsic surface density, the bins are 0.8 magni- 
tudes wide and are centered at 27.1, 27.9, and 28.7. The 
intrinsic apparent magnitude of each dropout was deter- 
mined by dividing the flux by the magnification at the 
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Figure 8. The source plane area effectively observed as a 
function of intrinsic Higo magnitude. This is the source 
plane area over which a galaxy with intrinsic magnitude mHi„t 
would be detected and satisfy our dropout selection criteria. 
The maximum possible area at a given intrinsic magnitude 
assuming 100% completeness is f2src = 1.57 arcmin'^. 



Errors in the number of objects in our surface den- 
sity and LF include both Poisson errors and sample vari- 
ance due to large scale structure, added in quadrature. 
The Poisson errors"'^'' are the 95% confidence intervals 
from the Bayesian method of Kraft et al. (1991). We 
follow Somerville et al. (2004) in calculating the sam- 

We use the sum of the fraction of the Icr error bars as the 
number of objects in each bin as discussed in Section 3. 
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pie variance and assume a power-law correlation func- 
tion ^ (r) ~ {ro/ry . Given the lack of constraints on 
the clustering of LBGs at high redshifts, we assume that 
the clustering of Lya emitters is sufficiently similar to 
inform our choice of 7 and rg (see e.g., Nagamine et al. 
2010). We fix 7 = 1.8 and assume tq = 7.0 h~^Mpc, 
which is conservative for z ^ 7 Lya galaxies according 
to Orsi et al. (2008, Figure 5). This gives a relative sam- 
ple standard deviation of ^ 38% for our 1.57 arcmin^ 
intrinsic solid angle, in agreement with the sample vari- 
ance calculator of Trenti & Stiavelli (2008) which gives 
37%. Given our counts, Poisson errors dominate over 
sample variance. 

The completeness corrected z ^ 7 counts of the 
HUDF09 blank field survey of Bouwens et al. (2011) are 
also shown in Figure 7. We expect our counts to be 
higher than those of Bouwens et al. (2011) as is the case. 
This is due to the fact that Bouwens et al. (2011) uses 
Z850! Y105 and J125 to search for zgso dropouts and given 
that Y105 and J125 combined cover the same wavelength 
range as our Jno, our redshift selection extends to higher 
redshift than that of Bouwens ct al. (2011). We do ob- 
tain counts to only slightly shallower depths than the 
Bouwens et al. (2011) deep blank field survey with dras- 
tically less observation time, demonstrating the benefit 
of cluster surveys in the search for high redshift galaxies. 
In addition, by searching for galaxies behind the Bul- 
let Cluster we see galaxies magnified by factors of 2-12 
(0.75-2.7 magnitudes). 

The completeness corrected LF of our dropout sample 
is shown in Figure 9 as a function of both lensed and 
intrinsic magnitudes. The best-fit z ^ 7 LF of Bouwens 
et al. (2011) is also shown. We first discuss the lensed LF. 
The bins are 0.8 magnitudes wide and centered at -21.5 
and -20.7. The lensed LF was calculated following Steidel 
et al. (1999). The lensed comoving volume effectively 
observed as a result of the completeness is a function of 
lensed Higo magnitude given by 

VeS,len (mHiao) = ^""3 J '^^'J^P ("^Hiso : z) : (H) 

where dz {dV/dz) is the comoving volume per unit solid 
angle between redshift z and z + dz and p (mHieo , z) is the 
probability that an LBG of lensed magnitude mHie,, 
redshift z will be detected in our images and satisfy our 
selection criteria (see Section 6). An estimate of the com- 
pleteness corrected comoving number density of galaxies 
per magnitude bin (the LF) is then given by the lensed 
surface density (before dividing by Aeffjen) divided by 
Vcffjon (niHieo) and the binsize. Vcff jen (mHiso ) was ap- 
plied to each bin by dividing each dropout within the 
bin by T4ffjon (niHieo) hnearly interpolated at the lensed 
magnitude of the dropout. The error bars include sam- 
ple variance and Poisson errors as described above. The 
fact that our lensed counts are significantly higher than 
the Bouwens et al. (2011) LF demonstrates the benefit of 
cluster surveys in this magnitude range as a result of the 
steepness of the LF leading to a positive magnification 
bias (see Section 1). 

We now discuss the intrinsic LF. The bins are 0.8 
magnitudes wide and centered at -19.8, -19, and -18.2. 
Kiff.int (niHieo int) the intrinsic comoving volume ef- 
fectively observed as a result of the completeness. It 
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Figure 9. Lensed (triangles) and intrinsic (circles) luminos- 
ity function of zgso dropouts with bins of width 0.8 magni- 
tudes. The lensed bins are centered at -21.5 and -20.7 and 
the intrinsic bins are centered at -19.8, -19, and -18.2. The 
Schechter parameterization of the best-fit luminosity function 
of the blank field study of Bouwens et al. (2011) is shown for 
comparison (solid line, cj)* = 0.9 x 10"^ Mpc"^, M* = -20.11, 
and a — —1.94). Uncertainty in our completeness estimate 
(see Section 6) is not reflected in the error bars. 

is a function of intrinsic magnitude calculated following 
Equation 11 with Jlimg replaced by rigic and p (muigoj z) 
replaced by p (mHigo int i z) . The intrinsic LF is then 
given by the intrinsic surface density (before dividing 
by Acff,i„t) divided by Feg^nt (mn 

160, int) and thfi binsize. 
Vcff.int (niHii3o,i„t) was apphed to each bin by applying it 
to each dropout within the bin given its intrinsic magni- 
tude. 

Our intrinsic LF is in agreement with that of Bouwens 
et al. (2011). It is not our intent to constrain the LF 
given so few objects, but we do provide further proof of 
the concept that cluster searches for high redshift LBGs 
are not only feasible but in many ways preferable. These 
cluster surveys also provide an independent check of high 
redshift blank field results at much lower observing cost. 
With a larger sample of massive clusters, with well under- 
stood magnification and magnification error properties 
(at least 2 multiply-imaged systems with known redshifts 
and weak-lensing data) and observed to similar depths 
as the Bullet Cluster, we can increase the number counts 
and significantly reduce the errors due to sample variance 
and Poisson sampling. 

8. CONCLUSIONS 

We find 10 candidate z ^ 7 zgso dropout galaxies 
behind the Bullet Cluster. Using the 8 objects with 
niHieo *^ 26.5, where the completeness is greater than 
22%, we calculate the surface density and luminosity 
function as a function of their intrinsic brightnesses. We 
find results consistent with published blank field surveys. 
We thus provide an independent check of blank field re- 
sults at only slightly shallower depths. The magnifica- 
tions of our Z850 dropouts range from 2-12 allowing us to 
detect sources up to 2.7 magnitudes deeper than blank 
field surveys with the same exposure time. With the 
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magnification of the Bullet Cluster we are thus able to 
probe to similar depths as blank field surveys despite us- 
ing much shallower data. 

In searching for high redshift LBGs, we find 8 (3) ob- 
jects with muieo < 25 that although they do not sat- 
isfy some of our zgso (Jno) dropout selection criteria and 
are thus rejected as low redshift contaminants, they are 
sufficiently red to satisfy our color cuts and are thus in- 
teresting in their own right. We give their positions and 
photometric properties. 

Magnifications are calculated from an optimally com- 
bined weak and strong lensing mass reconstruction of 
the Bullet Cluster (Bradac et al. 2009). Errors on the 
magnification are smaller than the Poisson sampling and 
sample variance. With more clusters with mass models 
of similar quality as the Bullet Cluster we could signifi- 
cantly reduce these errors and probe the luminosity func- 
tion at greater depths and higher redshifts than is other- 
wise possible. With a larger sample of efficient lenses we 
would also expect to find many highly magnified images, 
which would allow for better morphological studies and 
greatly ease the spectroscopic follow up of z > 7 galaxies. 
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